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Abstract 

We propose the Assembly Line Worker Integration and Balancing Problem (ALWIBP), a 
new assembly line balancing problem arising in lines with conventional and disabled workers. 
The goal of this problem is to maintain high productivity levels by minimizing the number 
of workstations needed to reach a given output, while integrating in the assembly line a 
number of disabled workers. Being able to efficiently manage a heterogeneous workforce is 
especially important in the current social context where companies are urged to integrate 
disabled workers. In this paper we present mathematical models and heuristics that can help 
assembly line managers to cope with this additional complexity; demonstrating by means 
of a robust benchmark how this integration can be done with losses of productivity that are 
much lower than expected. 

Keywords: Assembly line balancing; disabled workers; mathematical modeling. 



1 Introduction 



According to the International Labour Organization (ILO), people with disabilities represent 
an estimated 10 per cent of the world's population, including approximately 500 million of 
working age; being apparent that in the unemployment rates of the disabled are much higher 
than the average. 

Employment is the main path for social inclusion and participation in modern societies. 
Having a job is not only the basis for the survival and stability for many individuals, but 
also a key way of accessing many rights as citizens. Therefore the welfare and the social 
inclusion of the disabled depend very much on the degree of labor integration they are able to 
achieve. Different active policies to fight against discrimination have been set during the last 
decades, following models that are more/less inclusive depending on the local culture. Across 
specific national legislations, a general common formula is to reserve a share of workplaces 
in ordinary companies for people with disabilities. This share normally increases with the 
size of the company and, depending on the country legislation, usually goes from 2% to 5% 
of the jobs. 

Unfortunately, it is also a common phenomenon in many countries that this share is 
not always respected. Therefore it is clear that the solution should come not only by legal 
imposition, but mainly by overcoming the prejudices about the capabilities of the disabled, 
and by the genuine commitment of ordinary companies to include integration programs in 
their strategies. The aim of this paper is to contribute in making this commitment easier: 1) 
by providing the production managers with practical approaches that ease the integration of 
disabled workers in the production lines; 2) by demonstrating that, through the approaches 
proposed, the productivity of production systems suffers little (and often none) decrease. 

Once stated the great importance of integrating Disabled into the workforce of ordinary 
companies, we should make a brief introduction on some previous work inspired on the 
specific scenario of the so-called "Sheltered Work Centres for Disabled" (henceforth SWDs). 
SWDs are a special work formula legislated in many countries (with different variants) whose 
only difference from an ordinary company is that most of its workers must be disabled, 
and therefore they receive some institutional help in order to be able to compete in real 
markets. This labor integration formula has been successful in decreasing the former high 
unemployment rates of countries like Spain, and one of the strategies used by SWDs to 
facilitate the l abor integration has been the adoption of assembly lines. In this sense, 



Miralles et al.1 (120071 ) were the first to evidence how the integration of disabled workers 



in the productive systems can be done without losing, even gaining, productive efficiency 
through the use of assembly lines. This pioneer reference defined the so-called Assembly 
Line Worker Assignment and Balancing Problem (ALWABP) and demonstrated how the 
division of worker into single tasks becomes a powerful tool for making certain workers 
disabilities invisible. 



1.1 Contribution and outline of this work 

ALWABP was inspired in the SWDs reality, where the very high diversity of most of the 
workers and their limitations are the main characteristics. This scenario is quite different 
to the one of an ordinary company; where the aim is to efficiently integrate in the workforce 
just some workers to cope with the 2% to 5% of disabled workers legislation requirements. 
In this case the problem supposes much less diversity in the input data, and can also be 
stated with very different approaches with respect to the ALWABP, regarding the objective 
function, the hypothesis and model defined, and the kind of appropriate solution procedures. 
The aim of this paper is to introduce and analyze this new problem that has been named 
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"Assembly Line Worker Integration and Balancing Problem" (ALWIBP). Our study aims to 
answer specific requirements that normally arise in assembly lines of ordinary companies, 
where only few disabled workers have to be integrated, providing the production managers 
with practical tools that ease the integration of disabled workers in the most efficient manner. 
We propose new mathematical models for the problem as well as a constructive heuristic 
based on the similarities between the proposed problem and the so-called Simple Assembly 
Line Balancing Problem (SALBP). 

The remainder of this paper is structured as follows: in Section |2j we state a formal cod- 
ification of the new problem and some extensions, analyzing their practical implications and 
reviewing references of the literature with useful related approaches. Section [3] then presents 
the corresponding IP models for the proposed versions of the ALWIBP while Section [4] de- 
scribes a fast heuristic that has been developed to solve the problem. A experimental study 
in order to analyze the effectiveness of the proposed models and algorithms is conducted in 
Section [5j General conclusions end this manuscript. 



2 The Assembly Line Worker Integration and Balancing Prob- 
lem 

2.1 Introduction: SALBP vs ALWABP 

The SALBP was initially defined by Bavbarsl ( 1986 ) through several well-known simplifying 
hypotheses. This classical single-model problem, which consists of finding the best feasible 
assignment of tasks to stations so that certain precedence constraints are fulfilled, has been 
the reference problem in the literature in its two basic versions: when the cycle time C 
is given, and the aim is to optimize the number of necessary workstations, the problem is 
called SALBP-1. Whereas when there is a given number m of workstations, and the goal is 
to m inimize the cycle time C the literature knows this second version as SALBP-2 (jSchol] : 
1999h . 

A trend in Assembly Line research in the last decade has been to narrow the gap between 
the theoretical proposals and the industrial reality, which faces multiple specific configura- 



tions in detail such as multi-ma nned work s tations flD imitriadis; 2006; Ke llegoz and Toklu : 
2012 ). two sided assembly lines ( Kim et al. : 2009 : Ozcan and Toklul: 12009 ). or operator allo- 
cation in job sharing and operator revisiti ng lines ( Zeng et al. : 20121 ). among many others. 



As part of this trend, iMiralles et al.l (|2007l ) properly defined the ALWABP, a generalization 
of the SALBP where, in addition to the assignment of tasks to stations, a set of heteroge- 
neous workers also has to be assigned to stations. In this scenario each task has a worker- 
dependent processing time, which allows taking into account the limitations and specific 
production rates of each worker. Moreover, when the time to execute a task for certain 
worker is yery high, this assign ment is considered infeasible in the input data matrix. 

Since (IMiralles et al.l : 120071 ). many other references have contributed to give ALWABP 
visibility throughout academia, proposing different methods to solve the problem. The 
same authors have later developed a branch- and-bound algorithm for the problem, obtaining 



the exact solution of small-sized instances (jMiralles et al.l : 120081 ). Because of the problem 



complexity and the need to solve larger instances, the literature has since then shifted its 
efforts to heuristic methods. The current state-of-the -art methods for solving th e ALWABP 
are the iterated beam search (IB S) metaheurist i c of iBlum and Mirallesl (120111 ). the biased 
ra ndom-key ge n etic a lgorithm of iMoreira et al,l (|2012l ) and the iterative genetic algorithm 
of iMutlu et all ()2013l ). 
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2.2 ALWIBP definition 



The ALWABP problem was inspired in the SWDs reality with most workers presenting a 
high diversity of operation times; whereas the ALWIBP scenario introduced in Section 1 
pretends to simulate the "desirable" situation of having just a small percentage of disabled 
workers being integrated in a conventional assembly line. It has to be noted that the main 
(and only st udied) problem focusing on disabled i ntegr ation in assembly lines has been the 



ALWABP-2 (IMiralles et all ; 120081 ; iMoreira et all ; |2012| . e.g.), since the typical objective in 



SWD is to be as efficient as possible with the (diverse) available workforce. 

In the scenario associated with the ALWIBP, it makes sense to deal with the type 1 
problem, since the basic aim of a production manager can be to integrate the given disabled 
workers (in some cases some 2 or 5% of workers, or even more whether some compensation is 
needed due to low shares in other factory sections) while minimizing the number of additional 
workstations needed for doing so. This problem is named ALWIBP-1, by analogy with the 
SALBP case. 

In addition to this basic objective, once inside the solution subspace with minimal num- 
ber of workstations, the manager may aim to find those assignments in which the idle time 
in stations with disabled workers is minimum, in order to increase their partic ipation in the 



produ ction process. We call this extension ALWIBP-lS m i n . If, according to iBovsen et al. 



(120071 ) classification, ALWABP-2 was stated as [pa,link,cum\equip\c], in this case we can 
define ALWIBP-1 as [pa,link,cum\equip\m], while the ALWIBP-lS m i n can be stated as 
[pa, link, cum\equip\m, SLL stat ] using the same codification scheme. 

In the following, we propose integer linear models for the basic ALWIBP-1 situation and 
also for the extension proposed. 



3 Mathematical models 

In this section, we present a mathematical model for the ALWIBP-1 defined earlier, and 
further extend it to cope with ALWIBP 1-S m j n objectives with the use of the following 
notation: 



N set of tasks to be assigned; 

S set of workstations; 

W set of disabled workers, \W\ < \S\; 

ti execution time of task i when assigned to a conventional 

worker; 

t w i execution time of task i when assigned to disabled worker 

w G W; 

I-w Q N set of unfeasible tasks for worker w £ W; 
Fi set of immediate successors of task i; 

F* set of all successors of task i. 



3.1 Model for ALWIBP-1 

This f ormulation proposed for the ALWIBP-1 follows the idea used by Petterson and Albrachtl 
(| 19751 ) when modeling the SALBP-1. Let q be an artificial task and D q = {i£ N\F{ = 0} 
be the set of tasks that have no followers. We define a new precedence graph in which 
N' = N U {q} and all tasks in D q precede task q. The execution time of task q is always 0, 
i.e., t q = t wq = 0,Vu> € W. The ALWIBP-1 can thus be modelled as: 
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Mill ^ SX sq (1) 

ses 



subject to 



£a* = l, ViGiV', (2) 
^y sw = l, Vu> G W, (3) 

s£S 

Y Vsw < 1, Vs G 5, (4) 

^ x sl < ^ Vi, j € N',j eF i} keS,k^l, (5) 

s£S|s>fc s£S\s>k 

Y U ■ x si <C, Vs G 5, (6) 

^ • ^si < C 1 + L^(l - Vs G S,Vw G W, (7) 

ieN'\i w 

Usw < 1 - »st, Vs e5,VwG W, G (8) 

Y Vsw< Y x °n Vw€W,Vk€S,k^l, (9) 

seS\s>k seS\s>k 

x si e{0,i}, VseS,VieN', (10) 
y sw e{0,l}, VseS,Vw&W. (11) 

where: 

x s i binary variable equals to one if task i G N' is assigned to workstation 
s G 5, 

y su , binary variable equals to one if a disabled worker w G W is assigned to workstation 
a G 5, 

L w large constant, w G W. 

The objective function minimizes the index associated with the last station (the one that 
executes the fictitious last task q). In association with constraints (jSJ which state that all 
disabled workers are assigned, this objective function minimizes the number of conventional 
workers used in the line. Constraints garantee that each workstation receives at most 
one (disabled) worker. Constraints (|2J) ensure that all tasks are assigned, while constraints 
(jSJ) g uarantee that t he pre cedence relations are respected. These inequalities were proposed 



by iRitt and Costal ()2013f ) which analysed several versions of precedence constraints and 
concluded that constraints (jSJ presented better theoretical and practical results. 

Constraints (j5J) and (0 ensure that the cycle time is respected at stations without and 
with disabled workers, respectively. Constant L w must be sufficiently large to deactivate 
these last constraints if y sw = 0. Therefore, we take L w = YlieN\l w \^wi~ h\- This expression 
assumes the maximum additional time that a disabled worker w spends at a station in 
comparison to a conventional worker (this would be the additional time needed for the 
execution of all feasible tasks). 

Finally, constraints ([8j and (j5J) guarantee that disabled workers are not assigned to tasks 
which they are not able to execute and that they execute at least one task, respectively. 
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3.2 Extension for ALWIBP-1S 



In order to extend this formulation to the ALWIBP-lS m i n case, one can simply note that 
this new problem is characterized by the addition of another term in the objective function 
related to the idle time of the disabled workers. The new goal is to hierarchically minimize 
the number of stations (with higher priority) and the idle time of the stations with disabled 
workers. Thereby, this version of the problem aims to obtain more balanced solutions that 
increase the participation of these workers in production. 

To model this situation, we use non-negative real variables S w ,w G W, to measure the 
idle time of each disabled worker w. The values of these new variables are obtained with the 
aid of slack variables l sw , Vs G S, Vu; G W associated to constraints © , which are rewritten 
as: 

t wi - x si + l sw = C + L w (l - y sw ), Vs G S, Vu> G W, (12) 

i€N'\I w 

and with new constraints which are added to establish the correct relations between 5 W and 
the slack variables: 

6 W >lsw-(i-y sw )- [l w + ^ii) , VseS,VweW. (13) 



Variables l sw measure the slack for the constraints (jl 2 1) obtained with each pair s G 
S, w G W, while constraints (jl3[) define 5 W as the idle time of worker w, which is obtained 
with the value of l sw in the constraint corresponding to the actual assignment Usw 

The objective function can now include terms associated with these idle times, and the 
ALWIBP-lS m i n can be written as: 

Min Ysx sq + Y J^- (14) 

subject to 

© - ©,© - d, 

LeR+, Vs g S,yw G W, (15) 
^el + , Vs g S,Vw g W. (16) 

The constant term multiplying the idle time variables imposes a hierarchical character- 
istic in the objective function, giving priority to the minimization of stations and using the 
idle times as a secondary objective. 

4 Constructive Insertion Heuristic 

We propose a Constructive Insertion Heuristic (CIH) for the ALWIBP-1. This heuristic 
relies on the similarities between the ALWIBP and the classical SALBP. Indeed, during 
the heuristic procedures, a SALBP-1 solution is found and then iteratively adapted to 
incorporate the heterogeneous workers available in the ALWIBP-1. In the following section, 
this heuristic is presented in general terms, and then a formal description of the algorithm is 
presented in Section \4. 21 Finally, some variants of the procedure and two post-optimization 
routines are described in Sections 14.31 and 14. 4| respectively. 
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4.1 General structure of the proposed heuristic 

The main steps of the proposed heuristic are described in Algorithm [TJ 



Algorithm 1 Constructive insertion heuristic 
1: Ignore heterogeneous workers and obtain a SALBP-1 solution; 
2: Divide remaining line in segments (containing a certain number of stations each); 
3: Try all available heterogeneous workers in each station of the first line segment; 
4: Select best assignment and fix the solution on stations prior to the selected one. 
5: If there are still workers to assign, go to step 2; Otherwise, end. 



In the following subsections, each step of the algorithm is explained in details. 



4.1.1 Step 1: obtain a SALBP-1 solution 

A feasible solution for the SALBP-1, x re f, is used as starting point of the algorithm. This 
solution can be obtained with any heuristic approach. In this study, we use the best solution 
found by CPLEX 12.4 in 10 minutes of running time on a single thread of a 2.3Ghz Intel 
machine with 6Gb of memory. 



4.1.2 Step 2: divide line into segments 



The part of the line that has not yet been hxed is divided into \W a \ segments, where W a is 
the number of disabled workers not yet assigned. If a single worker is available, all stations 
in the remaining line are considered, otherwise, we consider stations in the current segment, 

m c - s b 



S r , defined as S r 



Sb 



l,...,s b + l + 



\W a \ 



where s& is the last fixed station and 



m r is the total number of stations in the current solution. 



4.1.3 Step 3: Try all possible assignments in current segment 

In this step, we try to obtain a solution by considering all assignments of workers in W a 
to stations in the current analyzed line segment S c , one at a time. When a station s € S c 
is considered for the assignment of worker w € W a , the following solution is obtained: in 
stations prior to s, the current solution, Xt m p, remains fixed. From s tation s and on, the solu- 



tion is rebuild using a constructive heuristic (jScholl and Vo/3l : Il996l ) . We call this procedure 
salbp-1 (x tmp ,w, s). The following heuristic criteria were used to pr ioritize the assignment 



of tasks during the procedure (the reader is referred to the article of IScholl and Vo/31 (| 19961 ) 
for more details): 

• MaxTime: Ascending order of task execution times, ti\ 

• MaxPW: Descending order of positional weights, •pw i = ij + YljzF? 

• MaxIF: Descending order of immediate followers, Fi] 

• MaxF: Descending order of followers, F*; 

In case a disabled worker is being considered, the appropriate parameters t w i and t w j 
are considered in criteria MaxTime and MaxPw. 
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4.1.4 Step 4: Select best assignment 



In this step, the best solution x* among those obtained for each assignment (w,s), with 
w € W a and s 6 S c is chosen. A solution x 1 obtained with an assignment of worker w\ to 
station s\ is considered better than a solution x 2 (x 1 < x 2 ) if the number of final stations 
of x 1 (designed by stations(x 1 )) is smaller than stations(x 2 ). If a tie happens, we consider 
the solution with larger idle time in the last station to be the best one. 



4.2 Formal structure of the proposed algorithm 

A formal description of the whole implemented heuristic is presented below. It uses a 
given SALBP-1 solution as input as described in Section [4.1.11 Line division into segments 
described in Section [4.1.21 is effected in lines 1-2 and in lines 20-21. The main loop of the 
algorithm occurs in lines 5-23: while there are workers to be assigned, each pair (w € W a , s £ 
S c ) is tested by obtaining a complete solution with the procedure described in Section 14.1.31 
(line 11) and the best one (according to the criterium presented in Section [4. 1.4|) is selected 
(lines 12-15). 



Algorithm 2 Constructive Insertion Heuristic 



Require: x re ^ : SALBP-1 reference solution; 
1: m c stations(x re *)\ 
m c 



2: S c 



3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17: 
18 
19 
20 

21: 



1,. 



Let x* be the incubent solution (initially, stations(x*) 
while W a ^ do 

Wb <— oo; 
S;, <— oo; 

for all w G W a do 
for all s£S c do 

Fix solutions of x re f on stations prior to s in x 
x c <— proc_salbp-l (x tmp ,w, s); 
if x c < x* then 
w w, 
Sb <- s; 
x* «- x c ; 
end if 
end for 
end for 

W a <- Wa\{w b }; 

m c <r— stations (x*); 



imp . 



re ^ <- X 



s b + 1, s b + 1 + 



Sb 



\Wa\ 



X 



22: 

23: end while 

Ensure: x*: best solution found. 



4.3 Algorithmic variants 



The proposed algorithm tries to assign workers from the first ordered station to the last 
one. A variant of this strategy might start at the last stations and proceed backwards. In 
this case, the line division into segments presented in Section 14.1.21 needs to be modified. 



Equations in line 2 and 21 of Algorithm [2] are modified for S c 



w a 



., m c > and 



7 



Sr. 



1 



\W a \ 



— 1 >, respectively. Since in this backward approach the tasks and 



workers assigned in the end of the line are kept fixed, the algorithm may need to insert 
intermediate stations. Therefore, the tie-breaker for the procedure described in Section [4.1. 41 
is no longer the idle time in the last station but the idle time in the last non-fixed station 
in the current solution. 

In addition to changing the order the ordered stations are visited for worker assign- 
ments, one can also change the order the tasks are assigned in the procedure described in 
Section 14.1.31 Indeed, as usual, the precedence graph can be considered in its normal or re- 
versed form, yielding four algorithmic variants (two worker assignment strategies combined 
with two task assignment strategies). Since the algorithm proposed is very fast, all four 
strategies were used at each run of the algorithm and the best result is kept; this increases 
the robustness of the algorithm with respect to the precedence graph structure. 



4.4 Post-optimization routines useful for ALWIBP-lS m j n 

The solution obtained by the CIH is improved by means of a mixed-integer programming 
neighborhood based on the original model (jl~])- (|lip . The main idea is to reduce the com- 
putational burden needed to solve the model by fixing the worker assignment variables and 
allowing changes only in task assignments variables. This can be done by solving the model 
with the addition of the following constraints: 

y s (w),w = 1, VwG W a , (17) 

where s(w) indicates the station to which worker w is assigned in the solution of the CIH. 
Clearly, model (fT])- (ITT]) can be adjusted accordingly, to eliminate useless constraints and 
variables. 

In order to further reduce the computational effort needed to solve the model, one can 
reduce the freedom of task assignment variables in the MIP neighborhood, by allowing tasks 
to be assigned only to the same station it was assigned in the CIH solution or to neighboring 
stations. This can be easily done with the addition of the following constraints: 



x a (i),i + aJ 8 (i)-i,i + x 9 (i)+l,i = h Vi G N, s(i) G S\{l,m}, (18) 
x u + x 2 i = l, V* G JV, s(i) = 1, (19) 
x m i + x m -i,i = 1, Vi G N,s(i) = m. (20) 

Where s(i) indicates the station to which task i was assigned in the CIH solution. 

Additionally, the objective function of this post-optimization program can be changed 
to incorporate the characteristics of the ALWABP-lS m i n , in which one prioritizes solutions 
with low idle times for the disabled workers. In this case, it suffices to use (|14p as objective 
function and add constraints: 

^2p s (w),i ■ Xsi + 5 W = C, Vu> G W, (21) 

Note that these simpler constraints replace constraints (|12p -(fT3|) in the original ALWIBP- 
ISmin problem, since the worker assignments are already known in the post-optimization 
phase. 
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5 Experimental study 



5.1 Justification of a new ALWIBP benchmark 

As discussed in section 2, the ALWABP was inspired in the situation found at SWDs where 
the very high diversity of workers and their limitations are the main characteristics; whereas 
the ALWIBP scenario pretends to simulate the "desirable" situation of only a small percent- 
age of disabled workers being integrated in conventional assembly lines. Moreover, as stated 
earlier, the main (and only studied) approach in this scenario has been ALWABP-2, since 
the typical objective in SWD is to be as efficient as possible with the (diverse) available 
workforce. Instead, in this new scenario, the assembly line balancing of type 1 (minimiza- 
tion of the number of stations given the desired cycle time) approach becomes realistic, 
since the basic aim of a production manager can be to integrate the normative (common 
in most countries) 2%-5% of disabled workers into the assembly line; while maintaining a 
given productivity. Many previous proposals fo r the ALWABP-2 w ere evaluated with the 
set of 320 benchmark instances first proposed by Chaves et al. ( 20091 ). Once stated the 



com- 



pletely different scenario where the ALWIBP arises, it is clear that this classical ALWABP 
benchmark is not useful here since, as explained above: (1) only a little share of the workers 
are disabled; and (2) the basic aim is now to minimize the number of workstations with 
non-disabled workers (ALWIBP-1 perspective). 

5.2 ALWIBP benchmark scheme 

As many other ALB approaches, the ALWABP benchmark was constructed from the only 
SALBP reference (the Scholl data collection of www.assembly-line-balancing.de); that was 
considered robust enough and has been extensively used to test mos t prop osals in the liter- 



ature so far. Nevertheless, as recently demonstrated bv lOtto et al.l (|2013l ). this framework 
does not seem rigorous enough. The problems were collected from different empirical and not 
empirical sources, and are based on only 25 precedence graphs; where just 18 distinct graphs 
have more than 2 5 tas ks and thus are meaningful for comparing solution methods. More- 
over, lOtto et al.l (J2013J) also point out the triviality of many classical instances, proposing 
to avoid these two inconvenient characteristics (which is: low diversity of graphs structure, 
and triviality), a SALBP generator and a new very robust challenging benchmark whose 
graphs morphologies include a sufficient variety of chains, bottlenecks and modules. Basi- 
cally, they propose different cells of data sets (with 25 different instances per cell) following 
a full-factorial design for the following parameters: 

1. Type of the graph: precedence graphs containing more chains, more bottlenecks, or 
mix of both. 

2. Order Strength: "low", "medium" and "high". 

3. Distribution of task times: "peak at the bottom", "bimodal" and "peak in the middle". 

4. Number of tasks: "small", "medium", "large" and "very large". 



Thus, we selected as basis the following collection of subsets from the lOtto et al.l (|2013f ) 
benchmark: 

1. We consider that diversity of graphs is sufficiently ensured selecting only the "mixed" 
instances (that have both chains and bottlenecks). 
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2. 



From them, we only select the data subsets with "low" and "high" Order Strength 
(that would give us clearer correlations if needed). 



3. And from them, regarding distribution of task times, we select the subsets with "peak 
at the bottom" and "bimodal" distribution (we discard the "peak in the middle" subset 
because the optimal number of stations is unknown for many of the instances). 

4. Finally, we take the resulting 100 robust instances of the "medium" (with n = 50 
tasks), "large" (with n = 100 tasks), and "very large" (with n = 1000 tasks) sub- 
sets, and we generate three corresponding benchmarks that we use separately in our 
experimentation (we discard the "small" subset since it is advised for simple testing 
only). 

In all three cases we generate the benchmark using the same procedure: for each 100 
instance we respect the precedence network and the conventional task time, and then we 
generate four different instances by adding one disabled worker with: high or low variability 
of task time respect to the original ones, and high or low percentage of incompatibilities. 
The two levels defined for the task times variability used the distributions U[ti, 2tj] and 
U[ti,5ti] for low and high variability, and the low and high percentage of incompatibilities 
in the tasks-workers matrix was set to 10% and 20% approximately. Following the same 
scheme we created 400 additional instances with two workers, then three workers, and finally 
four workers. 

Following this outline, we finally obtain three reliable benchmarks with the structure 
described, one with 1600 "medium" instances, a second one with 1600 "large" instances, 
and a third one with 1600 "very large" instances. It is important to note that in our source 
benchmark the instances are classified from "less tricky" to "extremely tricky" and it happens 
that, the four cells finally selected as base have a quite symmetric composition regarding 
the triviality. 

Finally, it is also important to note that the following ALWBIP experimental study has 
always used an in put cycle time of 1000 time units simply because it is the value used by 



Otto et al.l (|2013f ). and this way our results are comparable. 



5.3 ALWIBP experimental study 

Our experimental study consists of three parts: in the first experiment we use the "medium" 
and "large" benchmarks, obtaining exact solutions for almost all the instances with the model 
(and extension) of section 3. In the second experiment we use the same two benchmarks to 
compare the exact results with the ones obtained by the heuristic procedures. In the third 
experiment we use the "very large" benchmark to get further evidence of the good behavior 
of the heuristic against large problems that are difficult to solve exactly. 



5.3.1 Experiment 1: exact results for ALWIBP-1 and ALWIBP-lS min 

One basic aim of every company should be to integrate at least the normative percentage of 
disabled workers into the workforce. In this first experiment, we aim to demonstrate that the 
proposed methods enable the inclusion of higher percentages of disabled workers in the line 
without important losses in productivity. It has to be noted that productivity always means 
somehow (output result / input resources involved), and in assembly lines productivity can 
be defined as (throughput rate / number of workstations). Therefore (considering that in 
this experiment the Throughput rate is fixed), an increase on the number of workstations 
means, in general, a decrease of productivity. 
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Having this in mind, this first experiment was devoted to check what would be the shape 
of this expected loss of productivity. Both ALWIBP-1 and ALWIBP-lS m j n models could 
be solved in reasonable computation times using the commercial package CPLEX 12.4 (1 
thread, limit time of 1800 sec, and 6Gb as tree limit) with Intel Core i7 3.4 GHz, 16Gb 
RAM. 

In Tables [2] and [31 the columns indicate: 

• A: number of instances solved to optimality; 

• t(s): computational time (on average); 

• m-f and <7 mt : number of additional workstations needed in the ALWIBP-1 solution, 
with respect to the best know solution of SALBP-1 (average and deviation); 

• m-|- and o"m t (%) : percentage of additional workstations needed in the ALWIBP-1 
solution, with respect to the best know solution of SALBP-1 (average and deviation); 

• r: idle time of stations with disabled workers (on average) in ALWIBP-1; 

• r Smin : idle time of stations with disabled workers (on average) in ALWIBP-lS m j n ; 

• r)(%): percentage of tasks performed by disabled workers with respect to the mean 
number of tasks assigned to ordinary workers in ALWIBP-1; 

• r ?Smin(%) : percentage of tasks performed by disabled workers respect to the mean 
number of tasks assigned to ordinary workers in ALWIBP-lS m j n . 

• f3 and a /?(%): percentage of the number of disabled in the assembly line 
where m is the number of stations in the ALWIBP-1 solution (average and deviation); 

• 6: number of instances which there is no increse of the number of stations with respect 
to the SALBP-1 solution; no increase in the number of stations with respect to the 
SALBP-1 solution (average and deviation). 

We analyze first the results for "middle" instances. As expected, the increase in the 
number of stations grows with the number of disabled workers to be integrated and with 
the variability of the task times. Nevertheless, it can be observed that even in the most 
constrained case (4 disabled workers with execution times of up to 5 times the conventional 
time and 20% incompatibility), an average of only 1.6 new stations had to be added to 
integrate the workers. 

This is an even better result if we put the figures in context. In column 12, the f3 metric 
shows us that an average increase of 0.4 workstations imply in approximately 20% of the 
workforce composed by disabled workers. This is, in fact, relevant in real contexts, once the 
legislation of many countries requires a rate of disabled workers ranging from 2% to 5% in 
industries. 

Concerning the 6 criteria, we see that in 38% of the cases, we obtain a line balancing 
of ALWIBP-1 with the same number of workstations of the SALBP-1 solution, resulting in 
the best of possible scenarios (integrate all disabled workers without increasing production 
costs). Furthermore, the percentage of disabled workers that take part in the workforce 
keep at least in 5% in these cases. 

As the sub-space of possible ALWIBP-1 optimal solutions is large, we can combine 
this primary aim of minimizing conventional workstations with the (important) secondary 
objective of minimizing the idle time of disabled workers. After applying the extension 
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Table 2: Computational results: ALWIBP-1 and ALWiBP-lS min models (\N\ = 50) 



Var 



Inc 



t(s) I m-|- ± CTm t m-f- ± CT mt (%) 



»?(%) 



,(%) 



»±<t^(%) I 



U[t,2t] 
U[t,5t] 



10% 
20% 
10% 
20% 



95 104.0 

95 90.6 

95 125.0 

95 120.4 



0.2 ± 0.4 
0.2 ± 0.4 
0.4 ± 0.5 
0.4 ± 0.5 



2.8% ± 5.5% 
2.9% ± 5.6% 
5.0% ± 6.3% 
5.4% ± 6.5% 



149.9 
131.4 
231.6 
244.1 



1.6 
3.1 
9.6 
10.2 



66.2% 
62.8% 
43.7% 
39.9% 



87.6% 
85.6% 
63.4% 
58.9% 



10.9% 
10.9% 
10.6% 
10.6% 



± 2.7% 
± 2.7% 
± 2.6% 
± 2.5% 



78 
78 
59 
57 
63 
62 
30 
30 
48 
43 
8 
6 

32 
28 

1 





U[t,2t] 
U[t,5t] 



10% 
20% 
10% 
20% 



92 177.7 
96 85.8 

93 166.3 
93 152.7 



0.4 ± 0.5 
0.4 ± 0.5 
0.7 ± 0.5 
0.7 ± 0.5 



4.7% ± 6.5% 
4.7% ± 6.3% 
8.5% ± 6.8% 
8.6% ± 7.1% 



105.9 
101.3 
127.1 
140.1 



1.9 
3.5 
7.4 
8.1 



67.9% 
63.0% 
51.6% 
47.7% 



81.6% 
79.9% 
64.3% 
60.6% 



21.3% 
21.3% 
20.5% 
20.5% 



± 5.0% 
± 5.2% 
± 4.8% 
± 4.8% 



U[t,2t] 
U[t,5t] 



10% 
20% 
10% 
20% 



91 195.7 

89 213.6 
95 109.6 

90 217.0 



0.5 ± 0.5 

0.6 ± 0.5 

1.1 ± 0.5 

1.2 ± 0.5 



6.4% ± 6.5% 
6.9% ± 6.4% 
13.3% ± 8.1% 
14.1% ± 8.6% 



72.8 
76.5 
114.0 
126.5 



2.7 
3.1 
7.4 
11.0 



73.0% 
69.0% 
51.3% 
51.1% 



82.0% 
79.2% 
62.0% 
57.9% 



31.4% 
31.3% 
29.4% 
29.2% 



± 7.3% 
± 7.3% 
± 6.5% 
± 6.3% 



U[t,2t] 
U[t,5t] 



10% 
20% 
10% 
20% 



89 222.7 

86 301.2 

89 278.8 

84 362.9 



0.7 ± 0.5 

0.8 ± 0.5 

1.5 ± 0.6 

1.6 ± 0.6 



8.3% ± 6.8% 
9.2% ± 7.5% 
17.0% ± 9.7% 
18.1% ± 9.9% 



56.5 
68.0 
109.6 
120.5 



3.9 
5.0 
11.7 
11.7 



75.9% 
70.5% 
53.7% 
50.7% 



81.3% 
77.7% 
59.7% 
58.2% 



41.2% 
40.8% 
37.9% 
37.6% 



± 9.5% 
± 9.3% 
± 7.8% 
± 7.9% 



|W| 



Table 3: Computational results: ALWIBP-1 and ALWIBP-lS min models (\N\ = 100) 



Var 



Inc 



t(s) | m t ±<Tm T m 1 -±(Tm t (%) 



,(%) I P±<rp{%) I e 



U[t,2t] 
U[t,5t] 



10% 
20% 
10% 
20% 



86 
85 
81 
74 



309.1 
330.6 
435.2 
521.9 



0.2 ± 0.4 
0.2 ± 0.4 
0.3 ± 0.5 
0.4 ± 0.5 



1.0% 
1.0% 
2.0% 
2.4% 



± 2.3% 
± 2.3% 
± 2.9% 
± 3.1% 



68.9 
85.6 
109.1 
125.1 



1.1 
1.7 

3.4 
5.3 



74.4% 
67.0% 
56.8% 
49.7% 



87.0% 
85.7% 
66.0% 
67.0% 



5.6% ± 
5.6% ± 
5.5% ± 
5.5% ± 



1.4% 
1.4% 
1.3% 
1.3% 



83 
82 
66 
59 



U[t,2t] 
U[t,5t] 



10% 
20% 
10% 
20% 



67 675.8 

65 679.5 

53 971.1 

49 1045.3 



0.4 ± 0.5 

0.4 ± 0.5 

0.7 ± 0.5 

0.8 ± 0.4 



2.3% 
2.4% 
4.2% 
4.5% 



± 3.1% 
± 3.1% 
± 2.9% 
± 2.7% 



70.4 
72.7 
95.9 
105.5 



1.5 
1.7 
5.5 
4.8 



71.4% 
71.3% 
57.1% 
55.4% 



88.5% 
85.5% 
65.9% 
66.2% 



11.0% ± 
11.0% ± 
10.8% ± 
10.8% ± 



2.7% 
2.7% 
2.6% 
2.6% 



61 

59 

28 

21 
~W 

39 

5 

4 
"26" 

22 

2 

1 



U[t,2t] 
U[t,5t] 



10% 
20% 
10% 
20% 



47 1070.7 

47 1099.7 

28 1397.5 

28 1407.1 



0.6 ± 0.5 
0.6 ± 0.5 
1.1 ± 0.5 
1.1 ± 0.5 



3.3% 
3.5% 
6.4% 
6.5% 



± 3.1% 
± 3.0% 
± 3.2% 
± 3.2% 



48.8 
55.0 
86.7 
87.4 



1.7 
2.8 
5.5 
4.1 



74.6% 
73.4% 
57.4% 
58.0% 



85.2% 
83.8% 
65.1% 
63.3% 



16.4% ± 
16.4% ± 
15.9% ± 
15.9% ± 



3.9% 
3.9% 
3.7% 
3.7% 



U[t,2t] 
U[t,5t] 



10% 
20% 
10% 
20% 



34 1300.6 

32 1362.7 

8 1707.1 

12 1724.5 



0.8 ± 0.5 
0.8 ± 0.5 
1.5 ± 0.5 
1.5 ± 0.5 



4.4% 
4.6% 
8.4% 
8.5% 



± 2.9% 
± 2.9% 
± 4.0% 
± 3.7% 



55.2 
55.5 
89.3 
109.6 



2.7 
3.4 
7.4 
8.2 



74.1% 
72.4% 
57.7% 
57.6% 



82.6% 
80.8% 
61.4% 
58.5% 



21.6% ± 
21.6% ± 
20.8% ± 
20.8% ± 



5.1% 
5.1% 
4.7% 
4.8% 



ALWIBP-lS m j n to all the instances the results show that we maintain in all cases the same 
rough productivity (no increase of the number of workstations), while we reach a very little 
mean idle time (see column Ts min ). We note that in the column related to the rj(%) and 
^Smin(%) criteria, since the increase of these values follows the reduction of idle time in all 
scenarios studied. 

Concerning the ALWIBP models in the "large" instances (with 100 tasks), Table has 
shown that the obtention of optimal solutions becomes difficult when 3 and 4 disabled 
workers must be inserted in assembly lines (column A). However, the number of additional 
workstations is smaller than we notice in the experiments of "middle" instances, which 
implys that its approach may be useful in more complex assembly lines. Futhermore, the 
percentual of disabled workforce integrated continues in at least 5%, as it is desired. 

5.3.2 Experiment 2: Heuristics comparison 

For validation purposes, we solved the same 1600 "medium" and "large" instances with the 
CIH procedure and its variations with post-optimization phase (CIH-LS1 and CIH-LS2) 
presented in section 4. The MlP-local search routines were run with the package CPLEX 
12.4 with the same seeting previously shown, except the computational time limit, which 
was set to 60s. 

In Tables [4] and [5] the columns indicate: 

• mh t and cr mh : number of additional workstations in the CIH solution, with respect 
to the solution obtained when solving the ALWIBP-1 model (average and deviation); 

• mh t and cr mh (%): percentage of additional workstations in the CIH solution, with 
respect to the solution obtained when solving the ALWIBP-1 model (average and 
deviation) ; 

• th(s): Computational time used by the heuristic (on average); 

• Ties: Number of instances in which the CIH solution had the same number of work- 
stations (no increase) as the ALWIBP-1 model solution; 

• Ties sp i: Number of instances in which the CIH solution had the same number of 
workstations (no increase) as the SALBP-1 model solution; 

The results of Table |4]show a robust behaviour, where even in the worst case the increase 
of stations is 5,2% and the ties are 46. The computational times are practically null, which 
shows the computational efficiency of the heuristic. Concerning the CIH approaches with 
post-optimization routines, we see that the greater flexibility of the LSI causes significant 
improves in the results of CIH, increasing the number of ties compared with the solutions of 
ALWIBP-1 model in 24% (272 instances). We also outline that, in some cases, the average 
number of additional stations in the assembly lines is reduced by half. Despite the fact 
that LS2 improves the quality of solutions in only a few instances, it is very computationaly 
efficient, running in very little time (0.1s in average). 

Table [5] shows the results obtained for the large instances. We notice that the number 
of Ties (66%) is practically the same as observed in the previous case, which the average 
number of additional workstations is smaller. This states that the good results of this 
scenario is powered by the flexibility of assembly lines with many tasks. 

The computational times of the three approaches of the CIH show that the method is 
scalable. Moreover, the LSI increases the number of Ties related to the solutions of the 
ALWIBP-1 model, spending up to 53s on average, a very encouraging result if we consider 
that the CPLEX exceeds the computational time limite in most of cases. 
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Table 4: Computational results: CIH approaches (|iV| — 50) 



|W| 


Var 


Inc 




CIH 








CIH-LS1 








CIH-LS2 






m ht ± cr mht 


mh t ± o-m hT (%) 


t h (s) 


Ties 


m hf ± cr rah ^ 


m hT ± <r mhi . (%) 


t h (s) 


Ties 


m hf ± <7 rah ^ 


m hT ± <r rahi . (%) 


t h (s) 


Ties 




U[t,2t] 


10% 


0.2 ± 0.4 


1.8% ± 4.3% 


0.0 


84 


0.1 ± 0.2 


0.6% ± 2.7% 


10.5 


95 


0.1 ± 0.3 


1.3% ± 3.6% 


0.1 


88 


1 


20% 


0.1 ± 0.3 


1.3% ± 3.6% 


0.0 


88 


0.0 ± 0.1 


0.2% ± 1.7% 


10.8 


98 


0.1 ± 0.3 


1.1% ± 3.4% 


0.1 


89 




U[t,5t] 


10% 


0.2 ± 0.4 


2.1% ± 4.5% 


0.0 


81 


0.1 ± 0.3 


0.6% ± 2.6% 


10.8 


92 


0.2 ± 0.4 


1.8% ± 4.2% 


0.1 


83 




20% 


0.2 ± 0.4 


2.3% ± 4.7% 


0.0 


79 


0.0 ± 0.2 


0.4% ± 1.8% 


9.9 


96 


0.2 ± 0.4 


2.0% ± 4.4% 


0.0 


81 




U[t,2t] 


10% 


0.2 ± 0.4 


2.3% ± 4.6% 


0.0 


79 


0.1 ± 0.3 


0.9% ± 3.1% 


12.2 


91 


0.2 ± 0.4 


2.0% ± 4.4% 


0.1 


82 


2 


20% 


0.2 ± 0.4 


2.4% ± 4.7% 


0.0 


78 


0.1 ± 0.3 


1.0% ± 3.4% 


10.7 


91 


0.2 ± 0.4 


2.4% ± 4.7% 


0.1 


78 




U[t,5t] 


10% 


0.3 ± 0.5 


3.2% ± 5.1% 


0.0 


70 


0.2 ± 0.4 


1.7% ± 4.1% 


13.4 


85 


0.3 ± 0.5 


3.0% ± 5.0% 


0.1 


72 




20% 


0.3 ± 0.5 


3.4% ± 5.2% 


0.0 


69 


0.2 ± 0.4 


2.1% ± 4.4% 


12.1 


81 


0.3 ± 0.5 


3.4% ± 5.2% 


0.1 


69 




U[t,2t] 


10% 


0.3 ± 0.5 


3.3% ± 5.2% 


0.0 


69 


0.1 ± 0.3 


1.3% ± 3.4% 


16.4 


86 


0.3 ± 0.5 


3.1% ± 5.0% 


0.1 


71 


3 


20% 


0.3 ± 0.4 


2.8% ± 5.0% 


0.0 


74 


0.1 ± 0.3 


0.9% ± 3.0% 


14.1 


91 


0.2 ± 0.4 


2.4% ± 4.7% 


0.1 


78 




U[t,5t] 


10% 


0.4 ± 0.5 


3.9% ± 5.2% 


0.0 


62 


0.2 ± 0.4 


1.4% ± 3.5% 


17.0 


85 


0.4 ± 0.5 


3.9% ± 5.2% 


0.0 


62 




20% 


0.4 ± 0.5 


3.9% ± 5.1% 


0.0 


61 


0.2 ± 0.4 


1.6% ± 3.5% 


14.4 


83 


0.4 ± 0.5 


3.9% ± 5.1% 


0.0 


61 




U[t,2t] 


10% 


0.4 ± 0.5 


3.9% ± 5.4% 


0.0 


63 


0.2 ± 0.4 


2.0% ± 4.5% 


16.6 


82 


0.4 ± 0.5 


3.7% ± 5.3% 


0.1 


65 


4 


20% 


0.3 ± 0.5 


3.5% ± 5.1% 


0.0 


66 


0.1 ± 0.3 


1.3% ± 3.7% 


13.4 


84 


0.3 ± 0.5 


3.4% ± 5.1% 


0.0 


67 




U[t,5t] 


10% 
20% 


0.5 ± 0.5 
0.6 ± 0.5 


5.2% ± 5.1% 
5.2% ± 5.1% 


0.0 
0.0 


46 
46 


0.3 ± 0.4 
0.3 ± 0.5 


2.3% ± 4.2% 
2.6% ± 4.3% 


18.7 
14.6 


75 
72 


0.5 ± 0.5 
0.5 ± 0.5 


5.1% ± 5.1% 
5.2% ± 5.0% 


0.0 
0.0 


47 
46 











Table 5: 


Computational results: 


CIH approaches 


(W\ - 


100) 










|W| 


Var 


Inc 




CIH 








CIH-LSl 








CIH-LS2 






m ht ± cr mht 


m hf ± ^m hT (%) 


t h (s) 


Ties 


m hT ± <T rah ^ 


m hf ± (%) 


t h (s) 


Ties 


m hf ± <7 rah ^ 


m hT ± <T rah ^ (%) 


t h (s) 


Ties 




U[t,2t] 


10% 


0.3 ± 0.4 


1.5% ± 2.6% 


0.0 


73 


0.1 ± 0.3 


0.5% ± 1.6% 


47.5 


89 


0.2 ± 0.4 


1.1% ± 2.3% 


5.2 


79 


1 


20% 


0.3 ± 0.5 


1.8% ± 2.9% 


0.0 


71 


0.2 ± 0.4 


0.9% ± 2.3% 


46.9 


82 


0.2 ± 0.4 


1.3% ± 2.6% 


6.6 


77 




U[t,5t] 


10% 


0.4 ± 0.5 


1.9% ± 2.7% 


0.0 


64 


0.2 ± 0.4 


0.8% ± 2.0% 


43.4 


81 


0.4 ± 0.5 


1.9% ± 2.7% 


2.7 


65 




20% 


0.3 ± 0.4 


1.4% ± 2.5% 


0.0 


73 


0.1 ± 0.3 


0.4% ± 1.7% 


40.2 


87 


0.3 ± 0.4 


1.3% ± 2.4% 


2.8 


75 




U[t,2t] 


10% 


0.3 ± 0.4 


1.3% ± 2.4% 


0.0 


70 


0.2 ± 0.4 


0.8% ± 1.9% 


48.1 


80 


0.2 ± 0.4 


1.0% ± 2.2% 


2.3 


76 


2 


20% 


0.3 ± 0.4 


1.3% ± 2.4% 


0.0 


75 


0.2 ± 0.4 


0.8% ± 2.0% 


46.6 


82 


0.2 ± 0.4 


1.2% ± 2.3% 


2.1 


77 




U[t,5t] 


10% 


0.3 ± 0.5 


1.7% ± 2.7% 


0.1 


68 


0.2 ± 0.4 


0.9% ± 2.2% 


46.3 


74 


0.3 ± 0.5 


1.6% ± 2.6% 


3.0 


70 




20% 


0.3 ± 0.5 


1.6% ± 2.6% 


0.0 


72 


0.2 ± 0.4 


0.9% ± 2.2% 


43.9 


78 


0.3 ± 0.5 


1.6% ± 2.6% 


1.3 


72 




U[t,2t] 


10% 


0.3 ± 0.5 


1.5% ± 2.6% 


0.1 


70 


0.2 ± 0.4 


0.9% ± 2.1% 


49.6 


78 


0.3 ± 0.4 


1.4% ± 2.5% 


3.1 


72 


3 


20% 


0.3 ± 0.5 


1.6% ± 2.6% 


0.1 


71 


0.2 ± 0.4 


0.6% ± 1.9% 


46.9 


82 


0.3 ± 0.4 


1.4% ± 2.4% 


2.5 


74 




U[t,5t] 


10% 


0.4 ± 0.5 


2.1% ± 2.9% 


0.1 


63 


0.2 ± 0.4 


0.9% ± 2.1% 


53.4 


79 


0.4 ± 0.5 


2.1% ± 2.9% 


0.7 


63 




20% 


0.4 ± 0.5 


2.2% ± 2.8% 


0.1 


57 


0.2 ± 0.4 


0.9% ± 2.2% 


46.9 


75 


0.4 ± 0.5 


2.2% ± 2.8% 


1.2 


58 




U[t,2t] 


10% 


0.3 ± 0.5 


1.6% ± 2.5% 


0.1 


68 


0.3 ± 0.4 


1.0% ± 2.4% 


48.3 


68 


0.3 ± 0.5 


1.4% ± 2.4% 


1.7 


70 


4 


20% 


0.4 ± 0.5 


1.8% ± 2.8% 


0.1 


63 


0.2 ± 0.4 


0.9% ± 2.2% 


48.9 


74 


0.3 ± 0.5 


1.7% ± 2.7% 


2.0 


66 




U[t,5t] 


10% 
20% 


0.4 ± 0.5 
0.5 ± 0.5 


2.3% ± 2.8% 
2.7% ± 2.9% 


0.1 
0.1 


56 
52 


0.2 ± 0.4 
0.3 ± 0.4 


0.7% ± 2.2% 
1.3% ± 2.5% 


50.5 
49.6 


72 
70 


0.4 ± 0.5 
0.5 ± 0.5 


2.2% ± 2.8% 
2.5% ± 2.9% 


2.3 
0.9 


57 
55 



5.3.3 Experiment 3: ALWIBP experimentation with very large size instances 



The "very large" instances proposed by lOtto et al.l (120131 ) are much complex to be solved 
exactly by CPLEX. Although these instances are more theoretical than practical, compu- 
tational experiments with this scenario are important in order to put the scalability of the 
solving methods to proof. Thus, we highlight the ability of the CIH to solve problems of this 
dimension in a very reasonable computational times (3 min on average), as shown in Table 
[6] Additionally in 21% of the instances, the number of stations of a conventional assembly 
line did not change with the insertion of disabled workers. 



Table 6: Computational results: CIH (\N\ = 1000) 



|W| 


Var 


Inc 


CIH 


m ht ± a mht 


m ht ± CJm ht (%) 


t h (s) 


Ties S:p i 


1 


U[t,2t] 


10% 


0.8 ± 0.9 


0.5% ± 0.7% 


50.1 


34 






20% 


0.9 ± 0.9 


0.5% ± 0.7% 


50.8 


31 




U[t,5t] 


10% 


0.9 ± 0.9 


0.6% ± 0.7% 


50.4 


31 






20% 


1.0 ± 0.9 


0.6% ± 0.7% 


50.2 


29 


2 


U[t,2t] 


10% 


1.1 ± 0.9 


0.7% ± 0.7% 


106.6 


25 






20% 


1.1 ± 0.9 


0.7% ± 0.7% 


105.9 


24 




U[t,5t] 


10% 


1.5 ± 1.0 


1.0% ± 0.8% 


105.9 


16 






20% 


1.5 ± 1.1 


1.0% ± 0.8% 


105.9 


17 


3 


U[t,2t] 


10% 


1.2 ± 1.0 


0.8% ± 0.8% 


167.1 


21 






20% 


1.3 ± 1.0 


0.8% ± 0.7% 


167.3 


22 




U[t,5t] 


10% 


1.8 ± 1.1 


1.1% ± 0.8% 


164.4 


14 






20% 


1.9 ± 1.0 


1.2% ± 0.8% 


163.6 


8 


4 


U[t,2t] 


10% 


1.4 ± 1.0 


0.9% ± 0.8% 


223.4 


21 






20% 


1.4 ± 1.0 


0.9% ± 0.8% 


225.4 


21 




U[t,5t] 


10% 


2.1 ± 1.2 


1.4% ± 0.9% 


224.1 


10 






20% 


2.1 ± 1.2 


1.4% ± 0.9% 


222.0 


9 



6 Conclusions 

We propose the Assembly Line Worker Integration and Balancing Problem (ALWIBP), a 
new assembly line balancing problem arising in lines with conventional and disabled workers. 
This problem is relevant in a context where companies are urged to integrate disabled 
workers in their conventional productive schemes in order to cope with legislation issues or 
to include coorporate social responsibility goals in the production planning process. 

To solve this problem, we develop, first, an integer linear model that minimize the 
number of stations while ensuring the presence of all disabled workers in the assembly line. 
From this model, one variant that reduces the idle time of stations with disabled workers 
was proposed. We also implement a heuristic approach called CIH, that starts with a 
simple assembly line balancing situation and inserts the available disabled workers while 
reducing aditional stations. Furthermore, aiming the improvement of CIH solutions, two 
post-optimizations procedures based on MIP neighborhoods were designed. 

Through an experimental study on a extensive number of instances, we conclude that 
not only disabled workers can be included in the assembly lines with little productivity loss, 
but also that other planning goals can be simultaneously considered. Further work on this 
topic include the proposal of new adjacent objectives, implementation of more sophisticated 
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methods, and extensions that cope with job rotation schemes, U-shaped assembly lines and 
some lexocographic objectives. 
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